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Abstract 


LAMP  AT  and  ANSYS  are  analysis  tools  that  can  be  used  together  to  evaluate 
mechanical  integrity  of  structures  of  laminar  construction.  The  two  codes 
originate  from  two  different  areas  of  analysis.  Their  default  coordinate  systems 
have  different  orientations.  To  properly  examine  a  structure,  the  directional 
properties  used  by  LAMPAT  must  be  transformed  into  the  proper  input 
coordinate  directional  properties  for  ANSYS.  This  has  been  done  for  some  time 
already,  but  as  the  structures  being  analyzed  become  more  complex,  this 
coordinate  system  transformation  has  become  a  tool  requiring  some  additional 
thought.  This  report  describes  using  cylindrical  and  Cartesian  coordinate  system 
orientations  and  local  element  coordinate  system  results  within  ANSYS  to 
facilitate  the  analysis  of  laminate  structures  with  ply  layers  parallel  to  curvilinear 
profiles. 
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1.  Introduction 


To  properly  analyze  orthotropic  laminar  composites  in  structural  finite  element 
codes  such  as  ANSYS,  it  is  necessary  to  align  the  coordinate  systems  of  the 
laminae  with  the  coordinate  system  of  the  analysis  code.1  This  can  be 
accomplished  using  the  LAMP  AT2  code. 

Typically,  relatively  simple  laminate  geometries  have  been  analyzed  in  ANSYS, 
such  as  problems  with  structures  wound  on  cylindrical  mandrels  or  items 
machined  from  a  ply-stacked  composite  stock,  as  shown  in  Figure  1.  LAMP  AT 
uses  the  local  coordinate  system  reference  for  each  element  for  smearing  and 
unsmearing  the  element  data.  Models  of  these  simple  structures  can  be 
constructed  and  described  with  a  single  global  coordinate  system,  which  is  the 
local  coordinate  orientation  for  all  elements.  Thus,  it  is  relatively  easy  to 
properly  align  the  material  properties  for  analysis  in  ANSYS  and  LAMPAT.  The 
global  coordinate  system  may  be  cylindrical  or  Cartesian  (Figure  1).  ANSYS 
utilizes  the  element  coordinate  directions  and  the  smeared  properties  of  the 
elements  to  perform  its  analysis  on  an  approximated  laminate  structure. 
LAMPAT  analyzes  the  effect  of  the  stresses  on  each  element,  or  approximated 
laminate,  to  resolve  the  stresses  and  strains  in  the  individual  plies  within  that 
laminate.  The  material  properties  and  orientations  of  each  ply  are  contained  in 
the  LAMPAT  files. 

This  report  describes  how  to  use  LAMPAT  and  ANSYS  to  analyze  more  complex 
geometry  laminate  structures.  The  structure  in  Figure  2  represents  a  laminar 
structure  that  results  from  ply  stacking  into  a  mold  or  onto  a  pattern  of  complex 
geometry  to  produce  a  final  or  near  net  shape  part.  The  intralaminae  coordinate 
referencing  must  be  performed  as  it  was  previously,  but  additional  effort  is 
required  to  assure  that  the  laminate  is  properly  oriented,  and  the  resultant 
stresses  are  referenced  to  the  local  coordinate  system  for  each  element. 

Using  consistent  local  coordinate  directions  through  the  model  greatly  simplifies 
the  postprocessing  analysis  of  the  results.  Typically,  when  postprocessing  data 
from  LAMPAT  to  ANSYS,  the  failure  modes,  failure  criteria,  and  safety  factors  of 
the  various  regions  in  the  model  must  be  determined  to  properly  observe  the 
model's  response  to  the  loading.  For  example,  if  a  structure  like  the  one  in 


*  Swanson  Analysis  Systems,  Inc.,  P.O.  Box  65,  Johnson  Road,  Houston,  PA  15342-0065. 

1  Bogetti,  T.  A.,  C.  P.  R.  Hoppel,  and  B.  P.  Bums.  "A  Software  Tool  for  Analyzing  and  Designing 
Thick  Laminated  Composite  Structures."  ARL-TR-890,  U.S.  Army  Research  Laboratory,  Aberdeen 
Proving  Ground,  MD,  September  1995. 

2  Bogetti,  T.  A.,  C.  P.  R.  Hoppel,  and  W.  H.  Drysdale.  "Three-Dimensional  Effective  Property 
and  Strength  Prediction  of  Thick  Laminated  Composite  Media."  ARL-TR-911,  U.S.  Army  Research 
Laboratory,  Aberdeen  Proving  Ground,  MD,  October  1995. 


1 


GLOBAL  CARTESIAN  COORDINATE  SYSTEM 


Figure  1.  Typical  simple  laminate  structures. 


Figure  2.  A  laminate  structure  requiring  more  complex  coordinate  orientations  for 
analysis  purposes. 

Figure  3  is  being  analyzed  as  a  pressurized  container,  the  radial  and  hoop  loads 
of  the  structure  should  probably  be  observed.  If  the  conventional  approach  is 
used,  hoop  stress  will  be  SX  in  some  elements  and  SY  in  other  elements.  Radial 
stresses  will  behave  similarly.  Therefore,  it  will  not  be  sufficient  to  plot  SX  to  get 
radial  stress  and  SY  to  get  hoop  stress,  but  it  will  be  necessary  to  plot  the  failure 
mode  to  see  if  the  model  is  behaving  properly.  This  may  become  difficult  in 
structures  with  complex  shapes.  If  the  materials  are  oriented  so  that  the  Y  (and 
THETA  [0])  direction  for  all  of  the  elements  are  aligned,  plotting  SX  will  show 
radial  stress  for  all  elements,  and  SY  will  show  hoop  stress  for  all  elements. 

It  is  easier  to  align  the  coordinate  directions  of  the  materials  in  the  structure  and 
observe  the  structure's  behavior  in  a  few  plots,  than  it  is  to  keep  all  of  the 
regions'  results  and  conditions  and  put  them  together  like  a  puzzle  to  determine 
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GLOBAL  CARTESIAN  COORDINATE  SYSTEM 
IN  ANSYS 


Figure  3.  Alignment  of  local  coordinate  systems  in  ANSYS  to  facilitate  analysis  and 
postprocessing. 

the  structure's  behavior.  In  global  Cartesian  coordinates  in  ANSYS,  all  elements 
have  X,  Y,  and  Z  linear  directions;  stresses  in  postprocessing  will  be  relative  to 
this  X-Y-Z  order.  When  cylindrical  coordinates  are  used  either  globally  or 
locally,  the  order  will  be  R-0-Z.  Because  LAMP  AT  uses  the  element  coordinate 
systems  for  its  transformations,  it  may  be  important  to  define  them  to  be  other 
than  global  Cartesian  (the  default  system).  This  is  particularly  true  when  the 
global  Cartesian  system  cannot  properly  represent  the  system  being  analyzed. 
But  ANSYS  can  only  specify  an  element's  coordinate  system  as  global  Cartesian 
or  as  a  LOCAL  coordinate  system.  It  is  important  to  align  the  element- 
coordinate  systems  to  best  represent  the  structure  under  analysis.  Mixed 
coordinate  systems  should  be  aligned  to  simplify  postprocessing. 

This  report  describes  the  method  for  accomplishing  this  task.  Appendices  A-E 
include  example  input  files  for  implementing  this  task.  The  tools  required  for 
this  operation  include:  (1)  LAMP  AT  for  preprocessing  (smear)  the  material 
properties  for  analyzing  and  postprocessing  the  solution  (unsmear)  results  and 
writing  the  output  files,  (2)  ANSYS  for  solving  the  finite-element  problem,  and 
(3)  macros  (macro.write  and  macro.read)  for  writing  and  reading  files  from 
ANSYS  and  into  it. 
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2.  Method 


To  provide  a  simple  example  of  modeling  a  structure  with  complex  ply  contours, 
a  two-dimensional  (2-D)  cross  section  is  examined.  The  structures  can  be 
represented  as  a  series  of  arc  and  line  segments.  Several  local  cylindrical 
coordinate  systems  are  used  to  describe  these  arcs  in  ANSYS.  The  local 
cylindrical  coordinate  system  command  is  of  the  form: 

LOCAL/KCN,KCS,XC,YC,ZC,THXY,THYZ,THZX,PARl,PAR2/ 

where  KCN  is  the  reference  number  for  the  coordinate  system  (11-N), 

KCS  is  the  type  of  system  (0  Cartesian,  1  cylindrical,  2  spherical,  3  toroidal), 

XC,  XY,  and  XZ  are  the  global  Cartesian  coordinates  for  the  center  of  the  local 
system, 

THXY  is  the  rotation  about  local  Z  (X  towards  Y), 

THYZ  is  the  rotation  about  local  X  (Y  towards  Z), 

THZX  is  the  rotation  about  local  Y  (Z  towards  X),  and 
PARI  and  PAR2  are  used  for  describing  elliptical  shapes. 

The  direction  vectors  for  the  Cartesian  coordinate  systems  in  ANSYS  are  X,  Y, 
and  Z.  For  cylindrical  systems  they  are  R,  ©,  and  Z,  where  R  is  the  radial 
distance  from  the  center  of  the  system  in  the  YZ  plane,  0  is  the  angle  from  the  X 
axis  in  the  YZ  plane,  and  Z  is  the  distance  from  the  XY  plane. 

A  curved  cross  section  of  a  lamina  is  most  easily  described  with  a  cylindrical 
local  coordinate  system,  as  shown  in  Figure  4.  Thus,  the  out-of-plane  direction 
(3)  in  the  local  cylindrical  system  is  X.  The  zero  (1)  and  transverse  (2)  directions 
are  relative  to  0  and  Z,  which  depend  on  the  ply  orientation  within  the  laminae. 
For  a  unidirectional  ply,  the  1  direction  would  be  0,  the  2-direction  Z,  and  the 
3-direction  X  in  the  local  cylindrical  system.  For  a  continuous  layup  along  the 
pattern  surface  (Figure  3),  the  local  coordinate  systems  will  all  use  Y  or  0  as  the  1 
direction  of  the  ply.  This  is  because  the  descriptors  for  the  Cartesian  and 
cylindrical  coordinate  systems  are  given  in  the  orders  of  X,  Y,  Z  and  R,  0,  Z, 
respectively.  This  requires  that  local  coordinate  systems  be  constructed  as 
necessary  to  maintain  a  continuity  of  Y  and  ©  for  the  adjacent  elements  in  the 
two  reference  systems  (Figure  4). 

With  the  local  coordinate  systems  (LOCAL)  created  and  referenced  (CSYS)  to 
construct  the  proper  geometry  of  the  cross  section  and  continuity  of  the  plies, 
these  systems  must  be  referenced  using  the  ESYS  command  while  creating  the 
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Figure  4.  Curved  cross-section  laminate  structure. 

mesh  (i.e.,  using  an  AMESH  command).  ESYS,  KCN  assigns  a  coordinate 
reference  system  to  the  element.  The  default  state,  ESYS,  0,  assigns  global 
Cartesian  coordinates  to  the  elements.  The  user  may  input  (KCN)  either  0  for 
global  Cartesian  or  a  local  coordinate  system  (11-n).  The  local  coordinate  system 
must  have  been  previously  defined  using  the  LOCAL  command.  Neither  1, 2,  or 
3  (global  cylindrical,  spherical,  or  toroidal)  can  be  input  to  the  ESYS  command. 
Therefore,  even  for  a  model  that  can  be  built  or  described  using  the  global 
cylindrical  coordinate  system,  a  local  cylindrical  coordinate  system  with  its 
origin  at  0,  0,  0  must  be  defined  and  referenced  in  the  ESYS  command  to 
properly  define  the  material  properties  at  the  element  coordinate  system  level. 

Implementing  this  consistent  direction  reference  will  simplify  the  postprocessing 
of  the  results.  The  resultant  stresses  in  ANSYS  are  provided  by  default  and 
referenced  to  the  global  Cartesian  reference  system.  The  results  can  be  presented 
relative  to  other  coordinate's  systems  using  the  RSYS  command.  The  format  of 
the  RSYS  command  is  RSYS,  KCN.  KCN  can  be  0  for  global  Cartesian,  1  for 
global  cylindrical,  2  for  global  spherical,  11-n  for  local  coordinate  systems,  or 
SOLU  (in  menu  mode:  AS  CALCULATED)  for  reference  to  the  element 
coordinate  systems  defined  by  the  ESYS  command.  The  results  for  any  ANSYS 
analysis  can  be  viewed  (or  printed)  using  RSYS  (1  or  2,  or  3  or  11-n).  Only  by 
using  the  ESYS  command  can  the  RSYS,  SOLU  command  be  used  to  get  results  at 
the  same  time  relative  to  multiple  coordinate  systems.  If  the  ESYS  command  is 
not  implemented,  the  results  relative  to  the  single-selected  coordinate  system 
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will  be  applied  to  all  elements.  If  no  ESYS  command  is  applied,  the  elements 
have  their  default  elemental  coordinate  system,  which  is  global  Cartesian. 

The  following  steps  outline  a  procedure  for  conducting  a  laminate  structural 
analysis  with  ANSYS  and  LAMP  AT. 

(1)  Build  a  model  in  ANSYS  using  the  LOCAL,  CSYS,  and  ESYS  commands  to 
provide  consistent  fiber  orientation  within  each  ply  of  the  structure  (i.e., 
model.txt).  In  this  example,  the  model.txt  file  includes  all  prep7  input, 
except  for  the  material  property  data,  the  solution  section,  and  entry  into 
postl  including  selecting  the  results  set.  If  the  menu  mode  is  used,  the 
menu  paths  for  the  commands  may  be  obtained  from  the  help  menu.  Enter 
the  command  into  the  "Help  on"  window  in  older  versions  of  ANSYS  or 
enter  help,  "command  name,"  in  the  command  line  in  newer  versions  of 
ANSYS.  At  the  bottom  of  the  help  page  there  will  be  a  path(s)  to  the 
command  entry  point  in  the  menus. 

(2)  Create  a  material  input  deck  for  LAMP  AT  with  the  proper  ply  construction 
and  orientations  to  match  the  desired  directionality  in  the  ANSYS  model 
(i.e.,  1  =  Y,  2  =  Z,  3  =  X),  (matdb.in).  Note,  for  these  analyses,  that  the 
1  direction  is  considered  structurally  independent.  That  is,  for  example,  a 
hoop  winding  is  considered  a  "1"  direction  or  a  "0"  around  the  mandrel, 
not  a  90°  winding  as  in  conventional  processing  jargon. 

(3)  Rim  LAMPAT  as  the  preprocessor  to  create  a  material  property  input  deck 
from  matdb.in  to  use  in  ANSYS  (matiil). 

(4)  To  run  ANSYS  interactively,  start  ANSYS.  Click  on  Preprocessor.  Click  on 
Material  Props.  Click  on  Read  From  File.  Select  mat.fil.  Click  on  File. 
Click  on  Read  Input  From  ....  Select  model.txt.  In  this  example,  the 
model  will  be  run  and  solved,  and  ANSYS  will  be  waiting  in  the  postl 
postprocessor. 

(5)  To  get  data  for  postprocessing  in  LAMPAT  (you  are  still  in  ANSYS 
postprocessor),  click  on  File.  Click  on  Read  Input  From  ....  Select 
macro.write.  This  will  create  matidiil  and  stress.fil  output  files.  Copy  or 
move  these  files  where  they  can  be  used  with  LAMPAT. 

(6)  Using  the  postprocessor  mode  in  LAMPAT,  input  matdb.in,  matid.fil,  and 
stress.fil  when  prompted.  Use  lampat.fil  as  the  output  file.  (This  file  name 
will  be  used  later  to  read  data  into  ANSYS.)  Copy  or  move  lampat.fil  to 
where  it  can  be  accessed  in  ANSYS. 

(7)  In  ANSYS  in  postl,  click  on  File.  Click  on  Read  Input  From  ....  Select 
macro.read  (this  macro  reads  data  from  the  lampat.fil  file). 
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(8)  In  ANSYS,  execute  the  save  command  (either  enter  Save  in  the  command 
line,  or  select  the  SAVE  DB  button  with  the  mouse).  This  will  save  the 
etable  data  that  is  included  in  the  file.db  file,  not  in  the  result  file,  file.rst. 

This  process  creates  element  table  data  to  use  in  ANSYS  postprocessing.  The 
results  in  these  tables  will  be  relative  to  the  local  coordinate  system  for  each 
element.  Further  use  of  the  RSYS  command  in  postl  will  not  affect  the  element 
table  data  plotting  (Plot>Results>Contours>Elem  table),  but  it  will  affect  the  data 
plotted  in  the  ANSYS  results  files  (Plot>Results>Contours>nodal  or  elem). 

The  macro.write  and  macro.read  files  are  shown  in  Appendix  E.  They  were 
written  to  facilitate  the  data  transfer  between  ANSYS  and  LAMP  AT.  In  order  for 
all  of  the  LAMPAT  data  to  postprocess  in  ANSYS  via  etable  data,  both  the 
macro.write  and  macro.read  files  must  be  utilized.  The  resultant  etable  data  will 
be  SX,  SY,  SZ,  SYZ,  SXZ,  SXY,  SAFE-F,  FAIL-M,  and  CRIT-P.  Three  C4,  C5,  and 
C6  auxiliary  files  are  available. 


3.  Conclusions 


The  combined  use  of  LAMPAT  and  ANSYS  provides  a  very  powerful  tool  in 
analyzing  laminated  structures.  LAMPAT' s  use  of  element  coordinate  systems 
for  transforming  properties  and  results  requires  that  care  be  used  in  ANSYS  to 
properly  define  the  element  coordinate  systems  to  represent  the  structure  under 
analysis.  The  directional  properties  of  the  materials  must  be  properly  aligned  to 
the  structure.  The  local  coordinate  systems  in  ANSYS,  in  conjunction  with  the 
ANSYS  commands  of  ESYS  and  RSYS,  allows  for  the  proper  modeling  of 
laminate  structures  that  have  ply  layers  following  simple  or  complex  contours. 
With  proper  coordinate  system  rotations  and  ESYS  and  RSYS  commands, 
continuous  ply  orientation  along  curvilinear  profiles  can  be  more  easily  modeled 
and  postprocessed.  This  permits  the  stress  state  of  a  structure,  such  as  hoop  or 
radial  stress,  to  be  observed  in  a  single  plot  in  ANSYS  when  the  RSYS,  SOLU 
("RSYS  as  computed"  in  menu  mode)  command  has  been  issued.  ANSYS  cannot 
assign  the  global  cylindrical  as  an  element's  coordinate  system.  Even  when  the 
global  cylindrical  coordinate  system  can  be  utilized  to  construct  the  geometry,  a 
LOCAL  coordinate  cylindrical  system  (with  its  origin  at  0,  0,  0)  must  be  created 
to  use  as  input  to  the  ESYS  command.  The  ESYS  command  can  only  accept  input 
(KCN)  as  0  or  n,  where  n  is  an  integer  >10. 

When  analyzing  isotropic  materials  in  ANSYS,  the  ESYS  command  default  to 
global  Cartesian  is  sufficient  because  the  directional  properties  are  all  equal. 
Using  other  coordinate  system  results,  such  as  the  global  cylindrical  coordinate 
system  (RSYS,  1),  to  analyze  isotropic  materials  will  also  provide  the  proper 
results. 
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The  rotations  of  BETA  (B),  PHI  (O),  and  PSI  (¥)  correlate  to  the  rotation  of  Y 
towards  Z,  Z  towards  X,  and  X  towards  Y,  as  shown  in  Figures  A-l  and  A-2. 
LAMP  AT  solves  the  coordinate  rotations  only  in  the  order  of  BETA,  PHI,  and 
PSI,  so  there  are  a  limited  number  of  resultant  rotations.  The  permutations  are 
shown  in  Table  A-l.  The  values  of  2, 4,  and  7  are  provided  only  to  aid  in  tracking 
the  directional  changes  caused  by  the  rotations.  The  values  given  for  EX,  EY,  and 
EZ  represent  smeared  elastic  moduli  for  a  nonisotropic  material  calculated  by 
LAMP  AT.  The  BETA,  PHI,  PSI  (0,0,0)  state  reflects  the  results  of  the  calculated 
directional  values  based  on  the  percentage  of  fibers  from  the  various  fiber 
directions  inputted  for  each  laminate,  without  regard  to  the  laminate's 
orientation  in  the  structure  under  analysis.  The  laminate  is  oriented  to  the 
referenced  coordinate  system  with  the  B,  <£>,  and  'F  terms.  Users  might  configure 
a  simple  or  known  case,  such  as  a  unidirectional  lay  up,  and  make  a  trial  run  to 
understand  the  transformations  and  orientations  used  in  their  models.  It  may  be 
difficult  to  make  simple  checks  on  complex  geometry,  complex  lay  up  problems 
since  concepts  such  as  the  1,  2,  and  3  directions  become  muddled.  Figure  A-3 
shows  the  continuity  of  element  coordinate  direction  obtained  by  properly 
assigning  coordinate  systems. 


X  Y  Z 

REPRESENT  STANDARD 
COMPOSITE  CLAMPAT) 
COORDINATE  DIRECTIONS 

Xj  Y^Z 

REPRESENT  STANDARD 
ANSYS 

COORDINATE  DIRECTIONS 


X'  Y',Z' 

REPRESENT  THE 
LOCAL  ELEMENT 
COORDINATE  DIRECTIONS 


ROTATION  OF  BETA  E> 
MOVES  Y7  TOWARDS  Z7 


ROTATION  OF  PHI  <3> 
MOVES  Z7  TOWARDS  X' 


ROTATION  DF  PSI  ¥ 
MOVES  X'  TOWARDS  Y7 


Figure  A-l.  Coordinate  systems  in  LAMP  AT  and  ANSYS. 
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TO  GO  FROM  STANDARD  LAMPAT  TO  STANDARD  ANSYS  COORDINATE  SYSTEMS 
THE  ROTATION  IS  B  =  90,  <t>  =  90,  ¥  =  0 


TO  GO  FROM  STANDARD  LAMPAT  TD  CYLINDERICAL  ANSYS  COORDINATE  SYSTEMS 
THE  ROTATION  IS  B  =  90,  <f>  =  90,  ¥  =  0 


3  =  90  <p  =  90 


END  STATE 
IN  ANSYS 

9 


y?3  CSYS,0 
^  IN  ANSYS 


Figure  A-2.  Coordinate  transformations  for  consistency  within  ANSYS  and  LAMPAT. 


Table  A-l.  Permutations  of  coordinate  rotations. 


Orientation 

2 

3 

4 

5 

4 

3 

6 

BETA  (Y  to  Z) 

0 

90 

90 

90 

90 

0 

0 

0 

PHI  (Z  to  X) 

0 

0 

90 

90 

0 

90 

90 

0 

PSI  (X  to  Y) 

0 

0 

0 

90 

90 

0 

90 

90 

EX 

7 

7 

2 

2 

4 

2 

2 

7 

EY 

4 

2 

7 

4 

2 

4 

7 

4 

EZ 

2 

4 

4 

7 

7 

7 

4 

2 
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GLOBAL  CARTESIAN  COORDINATE  SYSTEM  GLOBAL  AND  LOCAL  COORDINATE  SYSTEMS 

ELEMENT  COORDINATE  SYSTEMS  ARE  NOT  ELEMENT  COORDINATE  SYSTEMS  ARE 

CONTINUOUS  AROUND  90  DEGREE  FILLET  CONTINUOUS  AROUND  90  DEGREE  FILLET 

Figure  A-3.  Comparison  of  default,  global  Cartesian,  and  local  coordinate  systems  on  the 
element  coordinate  systems  for  alignment  of  laminae  properties  and  results. 
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Intentionally  left  blank. 


Appendix  B.  Material  Input  for  Smearing  and  Coordinate 
Transformation  (DBMAT.IN) 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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Intentionally  left  blank. 
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8 

1 


NO.  OF  REGIONS  IN  THE  DATABASE 
NO.  OF  MATERIALS  IN  DATABASE 


1.0, 

2.0, 

o 

o 

H 

0.0 

!  REG_ID, 

FAIL_CRT, 

ITERS, 

ITEXT 

0.0, 

-90.0, 

0.0, 

0 . 0 

!  BETA_ (X) 

,  PHI__  ( Y)  , 

SI_(Z)  , 

— 

2.0, 

2.0, 

0.0, 

0.0 

!  NPLY,  REG_TYP,  - 

,  - 

1.0, 

0.005, 

0.0, 

-0 . 0001 

!  MAT_ID , 

THICKNESS, 

ANGLE, 

TEMP 

1.0, 

0.005, 

VO 

o 

o 

-0.0001 

!  MAT_ID, 

THICKNESS, 

ANGLE, 

TEMP 

2.0, 

2.0, 

o 

o 

H 

0.0 

!  REG_ID, 

FAIL_CRT, 

ITERS, 

ITEXT 

0.0, 

i 

vo 

o 

o 

0.0, 

0.0 

!  BETA_(X) 

,  PHI__(Y)  , 

SI__(Z)  , 

— 

2.0, 

2.0, 

0.0, 

0.0 

!  NPLY,  REG_TYP ,  --- 

,  - 

1.0, 

0.005, 

0.0, 

-0.0001 

!  MAT_ID, 

THICKNESS, 

ANGLE, 

TEMP 

1.0, 

0.005, 

90.0, 

-0.0001 

!  MAT__ID, 

THICKNESS, 

ANGLE, 

TEMP 

3.0, 

2.0, 

o 

o 

H 

0.0 

I  REG_ID, 

FAIL_CRT, 

ITERS, 

ITEXT 

0.0, 

-90.0, 

0.0, 

0.0 

!  BETA_ (X) 

,  PHI_  (Y)  , 

SI_(Z) , 

— 

2.0, 

2.0, 

0.0, 

0.0 

!  NPLY,  REG  TYP,  - 

,  - 

1.0, 

0.005, 

0.0, 

-0.0001 

!  MAT_ID, 

THICKNESS, 

ANGLE, 

TEMP 

1.0, 

0.005, 

90.0, 

-0.0001 

!  MAT_ID, 

THICKNESS, 

ANGLE, 

TEMP 

4.0, 

2.0, 

10.0, 

0.0 

!  REG_ID , 

FAIL_CRT , 

ITERS, 

ITEXT 

0.0, 

-90.0, 

0.0, 

0.0 

!  BETA  (X) 

,  PHI  (Y)  , 

SIJZ)  , 

— 

2.0, 

2.0, 

0.0, 

0.0 

I  NPLY,  REG_TYP,  --- 

,  - 

1.0, 

0.005, 

0.0, 

-0.0001 

!  MAT_ID, 

THICKNESS, 

ANGLE, 

TEMP 

1.0, 

0.005, 

VO 

o 

o 

-0.0001 

!  MAT_ID, 

THICKNESS, 

ANGLE, 

TEMP 

5.0, 

2.0, 

10 .0, 

o 

o 

!  REG_ID, 

FAIL_CRT, 

ITERS, 

ITEXT 

o 

o 

VO 

o 

o 

o 

o 

cr\ 

0 . 0 

!  BETA__  (X)  , 

PHI_JY)  , 

SI_(Z)  , 

— 

2.0, 

2.0, 

0.0, 

0.0 

!  NPLY, 

REGJTYP, 

— 

1.0, 

0.005, 

0.0, 

-0.0001 

!  MAT_ID , 

THICKNESS, 

ANGLE, 

TEMP 

1.0, 

0.005, 

90.0, 

-0.0001 

!  MAT_ID, 

THICKNESS, 

ANGLE, 

TEMP 

6.0, 

2.0, 

10.0, 

0.0 

!  REG_ID, 

FAIL_CRT, 

ITERS, 

ITEXT 

VO 

o 

o 

90.0, 

o 

o 

o\ 

0.0 

!  BETA_(X), 

PHI_  ( Y)  , 

SI_(Z)  , 

— 

2.0, 

2.0, 

0.0, 

0.0 

!  NPLY, 

REGJTYP, 

— 

o 

H 

0.005, 

0.0, 

-0.0001 

I  MAT_ID, 

THICKNESS, 

ANGLE, 

TEMP 

1.0, 

0.005, 

90.0, 

-0.0001 

»  MAT_ID, 

THICKNESS, 

ANGLE, 

TEMP 

1.0, 

2.0, 

1 . 889E-04 

t 

MAT_ID, 

MAT_TYPE , 

MAT  DENS (S2- GLASS/ 

7 . 1500E+06 , 

2.13E+06, 

2.13E+06 

t 

El, 

E2 , 

E3  EPOXY) 

0.499, 

0.306, 

0.296 

i 

NU23  , 

NU13  , 

NU12 

0 .710E+06, 

0 .980E+06, 

0 . 980E+06 

i 

G23 , 

G13, 

G12 

2 . 300E-06 , 

1 . 850E-05 , 

1.850E-05 

t 

Al, 

A2 , 

A3 

0.0, 

o 

o 

o 

o 

!  (1) 

VM1, 

VM2, 

VM3  VON  MISES 

0.0, 

0.0, 

0.0 

j 

VM4, 

VM5, 

VM6 

0.0, 

0.0, 

0.0 

t 

VM7 , 

VM8 , 

VM9 

o 

o 

0.0, 

0.0 

T 

- , 

- , 

— 

0.0, 

0.0, 

0.0 

t 

— 

243 . O0E+03 , 

7.00E+03, 

8.50E+03 

!  (2) 

X1T, 

X2T, 

X3T  MAXIMUM 

177 . O0E+03 , 

3  0 . 60E+03 , 

35 . 0OE+O3 

j 

X1C, 

X2C, 

X3C  STRESS 

15 . 70E+03 , 

17 . 0OE+03 , 

15 . 70E+03 

i 

X23, 

XI 3 , 

XI 2 
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0.0, 

0.0, 

0.0 

1 

- , 

— 

o 

o 

0.0, 

0.0 

1 

- / 

— 

0.0340, 

0.0033, 

0.0040 

!  (3) 

Y1T, 

Y2T, 

Y3T  MAXIMUM 

0.0248, 

0 . 0144  , 

0.0164 

i 

Y1C, 

Y2C, 

Y3C  STRAIN 

0 . 0221, 

0.0173, 

0.0160 

i 

Y23 , 

Y13 , 

Y12 

0.0, 

0.0, 

0.0 

t 

— 

o 

o 

o 

o 

o 

o 

i 

— 

0.0, 

0.0, 

0.0 

!  (4) 

XLT, 

XTT, 

HPD  HYDROSTATIC 

0.0, 

0.0, 

0.0 

j 

XLC, 

XTC, 

S  PRESSURE 

0 .0, 

o 

o 

0 . 0 

t 

ML1 , 

ML2 , 

LTP 

0.0, 

0.0, 

0.0 

i 

MT1 , 

MT2 , 

TTP 

0.0, 

0.0, 

0.0 

t 

MSI, 

MS2 , 

STP 

o 

o 

0.0, 

o 

o 

!  (5) 

X1T, 

X1C, 

-  TSAI ~ WU 

0.0, 

0.0, 

0.0 

i 

X2T, 

X2C, 

— 

0.0, 

0.0, 

o 

o 

i 

S12 , 

FI  2 , 

F23 

0.0, 

0.0, 

0.0 

i 

— 

0.0, 

0.0, 

0.0 

i 

— 

1500E+06 , 

2.13E+06, 

2 . 13E+06 

!  (6) 

El, 

E2 , 

E3  CHRISTENSEN 

0.499, 

0.306, 

0.296 

i 

NU23, 

NU13 , 

NU12 

.  710E+06 , 

0 . 980E+06 , 

0 . 980E+06 

j 

G23 , 

G13 , 

G12 

o 

o 

o 

o 

o 

o 

? 

K, 

ALPHA, 

— 

0.0, 

0.0, 

0.0 

t 

Y1T, 

Y1C, 

— 

o 

o 

o 

o 

o 

o 

!  (7) 

FE1 , 

FE2 , 

FE3  FENG 

o 

o 

0.0, 

o 

o 

j 

FE4  , 

FES, 

FE6 

0.0, 

0.0, 

0.0 

j 

FE7 , 

FE8 , 

FE9 

o 

o 

0.0, 

o 

o 

i 

FE10 , 

FE11 , 

FE12 

0.0, 

0.0, 

0.0 

i 

FE13 , 

FE14  , 

FE15 

0.0, 

0.0, 

0.0 

I  (8) 

HA1 , 

HA2 , 

HA3  MODIFIED 

0.0, 

o 

o 

o 

o 

i 

HA4  , 

HAS, 

HA 6  HASHIN 

0.0, 

0.0, 

0.0 

! 

HA  7 , 

HA  8  , 

HA  9 

0.0, 

0.0, 

o 

o 

J 

HA10 , 

HA11 , 

HA12 

o 

o 

0.0, 

o 

o 

1 

HA13, 

HA14, 

HA15 
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Appendix  C.  Material  Input  for  ANSYS  (MAT.FIL) 


This  appendix  appears  in  its  original  form,,  without  editorial  change. 
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Intentionally  left  blank. 


/COM,  LAMPAT  (VERSION  4.1)  MATERIAL  FILE  FOR  ANSYS  5.5 
MPTEMP 

MPTEMP,  1,  0 . 0000000E+00, 

MPDATA , EX  ,  1,1,  0 . 1130000E+07 , 

MPTEMP 

MPTEMP,  1,  0 . 0000000E+00 , 

MPDATA, EY  ,  1  ,  1,  0 . 7150000E+07 , 

MPTEMP 

MPTEMP,  1,  0 . 0000000E+00 , 

MPDATA, EZ  ,  1,1,  0 . 2130000E+07 , 

MPTEMP 

MPTEMP,  1,  0 . 0000000E+00 , 

MPDATA, NUXY,  1  ,  1,  0 . 3060000E+00 , 

MPTEMP 

MPTEMP,  1,  0 . 0000000E+00 , 

MPDATA, GXY  ,  1,1,  0 . 9800000E+06 , 

MPTEMP 

MPTEMP,  1,  0 . 0000000E+00 , 

MPDATA, GYZ,  1  ,  1,  0 . 9800 00 0E+06 , 

MPTEMP 

MPTEMP,  1,  0 . 0000000E+00 , 

MPDATA, GXZ  ,  1,1,  0 . 7099999E+06 , 

MPTEMP 

MPTEMP,  1,  0 . 0000000E+00 , 

MPDATA, PRXY,  1  ,  1,  0 . 4836084E- 01 , 

MPTEMP 

MPTEMP,  1,  0 . OOOOOOOE+OO, 

MPDATA, PRYZ,  1  ,  1,  0 . 2960001E+00 , 

MPTEMP 

MPTEMP,  1,  0 . 0000000E+00, 

MPDATA, PRXZ,  1  ,  1,  0 . 2647277E+00 , 

MPTEMP 

MPTEMP,  1,  0 . 0000000E+00, 

MPDATA, EX  ,  2,1,  0 . 1130000E+07 , 

MPTEMP 

MPTEMP,  1,  0 . 0000000E+00, 

MPDATA, EY  ,  2,1,  0 . 7150 00 0E+07 , 

MPTEMP 

MPTEMP,  1,  0. 0000000E+00, 

MPDATA, EZ  ,  2,1,  0 . 2130000E+07 , 

MPTEMP 

MPTEMP,  1,  0 . 0000000E+00, 

MPDATA, NUXY,  2  ,  1,  0 . 3060000E+00 , 

MPTEMP 

MPTEMP,  1,  0 . OOOOOOOE+OO, 

MPDATA, GXY  ,  2,1,  0 . 9800000E+06, 

MPTEMP 

MPTEMP,  1,  0 . 0000000E+00, 

MPDATA, GYZ,  2  ,  1,  0 . 9800000E+06 , 

MPTEMP 

MPTEMP,  1,  0 . 0000000E+00, 

MPDATA, GXZ  ,  2,1,  0 . 7099999E+0 6 , 

MPTEMP 

MPTEMP,  1,  0 . 0000000E+00, 


MPDATA,  PRXY, 


MPTEMP 

MPTEMP ,  1, 

MPDATA,  PRYZ , 
MPTEMP 

MPTEMP,  1, 

MPDATA,  PRXZ, 
MPTEMP 

MPTEMP,  1, 

MPDATA,  EX  , 
MPTEMP 

MPTEMP,  1, 

MPDATA,  EY  , 
MPTEMP 

MPTEMP,  1, 

MPDATA, EZ  , 
MPTEMP 

MPTEMP ,  1 , 

MPDATA,  NUXY, 
MPTEMP 

MPTEMP,  1, 

MPDATA, GXY  , 
MPTEMP 

MPTEMP,  1, 

MPDATA, GYZ, 
MPTEMP 

MPTEMP,  1, 

MPDATA, GXZ  , 
MPTEMP 

MPTEMP,  1, 

MPDATA, PRXY, 
MPTEMP 

MPTEMP,  1, 

MPDATA,  PRYZ, 
MPTEMP 

MPTEMP,  1, 

MPDATA, PRXZ , 
MPTEMP 

MPTEMP ,  1 , 

MPDATA, EX  , 
MPTEMP 

MPTEMP,  1, 

MPDATA, EY  , 
MPTEMP 

MPTEMP,  1, 

MPDATA, EZ  , 
MPTEMP 

MPTEMP,  1, 

MPDATA , NUXY , 
MPTEMP 

MPTEMP,  1, 

MPDATA, GXY  , 
MPTEMP 

MPTEMP,  1, 


2  ,  1,  0 -4836084E-01, 

0 . OOOOOOOE+OO, 

2  ,  1,  0 . 2960001E+00 , 

0 . 00000  0  0E+0  0  , 

2  ,  1,  0.2647277E+00, 

O.OOOOOOOE+OO, 

3  ,  1,  0 . 113  000  0E+07 , 

O.OOOOOOOE+OO, 

3  ,  1,  0 . 7 150  00  0E+07 , 

O.OOOOOOOE+OO, 

3  ,  1,  0 . 2 13  00  00E+07 , 

0 . OOOOOOOE+OO , 

3  ,  1,  0 . 3  06000  0E+00 , 

0 . OOOOOOOE+OO, 

3  ,  1,  0 . 9800000E+06, 

0 . 0000000E+00, 

3  ,  1,  0 . 9800000E+06 , 

0 . OOOOOOOE+OO, 

3  ,  1,  0 . 7099999E+06 , 

0.0000000E+00, 

3  ,  1,  0 .4836084E-01, 

0 . OOOOOOOE+OO, 

3  ,  1,  0 . 2  960001E+00 , 

0.0000000E+00, 

3  ,  1,  0 . 264  72 7 7 E+ 00  , 

0 . 0000000E+00, 

4  ,  1,  0 . 113  0  00  0E+07 , 

0  .  OOOOOOOE+OO, 

4  ,  1,  0 . 715000  0E+0 7 , 

O.OOOOOOOE+OO, 

4  ,  1,  0 .2130000E+07, 

0.0000000E+00, 

4  ,  1,  0.3  0600  0  0E+0  0 , 

0 . OOOOOOOE+OO, 

4  ,  1,  0 . 9800000E+06 , 

0. OOOOOOOE+OO, 
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4 


1,  0 . 9800000E+06 , 


mpdata,  gyz, 

MPTEMP 

MPTEMP  ,  1, 

MPDATA,  GXZ  # 
MPTEMP 

MPTEMP,  1, 

MPDATA ,  PRXY, 
MPTEMP 

MPTEMP ,  1 , 

MPDATA, PRYZ , 
MPTEMP 

MPTEMP,  1, 

MPDATA, PRXZ, 
MPTEMP 

MPTEMP,  1, 

MPDATA, EX  , 
MPTEMP 

MPTEMP,  1, 

MPDATA, EY  , 
MPTEMP 

MPTEMP,  1, 

MPDATA, EZ  , 
MPTEMP 

MPTEMP,  1, 

MPDATA ,  NTJXY , 
MPTEMP 

MPTEMP,  1, 

MPDATA, GXY  , 
MPTEMP 

MPTEMP,  1, 

MPDATA, GYZ, 
MPTEMP 

MPTEMP,  1, 

MPDATA, GXZ  , 
MPTEMP 

MPTEMP,  1, 

MPDATA, PRXY, 
MPTEMP 

MPTEMP,  1, 

MPDATA, PRYZ, 
MPTEMP 

MPTEMP ,  1 , 

MPDATA, PRXZ , 
MPTEMP 

MPTEMP,  1, 

MPDATA, EX  , 
MPTEMP 

MPTEMP,  1, 

MPDATA, EY  , 
MPTEMP 

MPTEMP,  1, 

MPDATA, EZ  , 
MPTEMP 

MPTEMP,  1, 


O.OOOOOOOE+OO, 

4  ,  1,  0 . 7099999E+06 , 

0. OOOOOOOE+OO, 

4  ,  1,  0 . 4  83  6084E-01 , 

0 . OOOOOOOE+OO, 

4  ,  1,  0 . 2960001E+00 , 

O.OOOOOOOE+OO, 

4  ,  1,  0 . 2  647277E+00 , 

O.OOOOOOOE+OO, 

5  ,  1,  0. 1130000E+07, 

0 . OOOOOOOE+OO, 

5  ,  1,  0 . 7150  000E+07 , 

0 .OOOOOOOE+OO, 

5  ,  1,  0 . 2130000E+07 , 

O.OOOOOOOE+OO, 

5  ,  1,  0 .3060000E+00, 

O.OOOOOOOE+OO, 

5  ,  1,  0 . 9800000E+06 , 

O.OOOOOOOE+OO, 

5  ,  1,  0 . 9800000E+06, 

O.OOOOOOOE+OO, 

5  ,  1,  0 . 7099999E+06 , 

O.OOOOOOOE+OO, 

5  ,  1,  0 . 4836084E-01 , 

O.OOOOOOOE+OO, 

5  ,  1,  0 .2960001E+00, 

O.OOOOOOOE+OO, 

5  ,  1,  0 . 2  64  7277E+00 , 

O.OOOOOOOE+OO, 

6  ,  1,  0.1130000E+07, 

O.OOOOOOOE+OO, 

6  ,  1,  0 .7150000E+07, 

O.OOOOOOOE+OO, 

6  ,  1,  0 . 2130000E+07 , 

O.OOOOOOOE+OO, 


MPDATA,  NUXY, 
MPTEMP 

MPTEMP ,  1, 

M  PD  AT  A,  GXY  , 
MPTEMP 

MPTEMP ,  1 , 

M  PD  AT  A,  GYZ  , 
MPTEMP 

MPTEMP,  1, 

M  PD  AT  A,  GXZ  , 
MPTEMP 

MPTEMP,  1, 

MPDATA , PRXY , 
MPTEMP 

MPTEMP,  1, 

MPDATA, PRYZ, 
MPTEMP 

MPTEMP,  1, 

MPDATA, PRXZ, 


6  ,  1,  0.3060000E+00, 

0 . OOOOOOOE+OO , 

6  ,  1,  0 . 9800000E+06, 

0 . OOOOOOOE+OO, 

6  ,  1,  0 . 9800000E+06 , 

0 . OOOOOOOE+OO, 

6  ,  1,  0.7099999E+06, 

0 .OOOOOOOE+OO, 

6  ,  1,  0.4836084E-01, 

0. OOOOOOOE+OO, 

6  ,  1,  0 . 2  960  0  01E+0  0 , 

0 . OOOOOOOE+OO, 

6  ,  1,  0.2647277E+00, 
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Appendix  D.  ANSYS  Input  File  (MODEL.TXT) 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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Intentionally  left  blank. 
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If  using  the  menu  mode,  the  menu  paths  for  the  commands  may  be 
obtained  through  the  help  menu.  Enter  the  command  into  the  "Help  on" 
window;  a  path(s)  to  the  command  entry  point  in  the  menus  will  be  at  the 
bottom  of  the  help  page. 


/title,  LAYUP:  0/90 
!  DEFINE  MATERIAL  PROPERTIES 

!  CHOOSE  ELEMENT  TYPE  2-D  STRUCTURAL 

et,l,42, , , 0 

et,2,42, , ,0 

et , 3 , 42 , ,  ,  0 

et , 4 , 42 , , ,0 

et,5,42, , ,0 

et,6,42, , ,0 


!  CREATE  LOCAL  COORDINATE  SYSTEMS 

I  MENU  PATH  FOR  LOCAL:  Utility  Menu  >WorkPlane  >Local  Coordinate  Systems 
!  >Create  Local  CS 


local, 11, 1,-1, -1,0 
local, 12, 1,1, -1,0 
local, 13, 1,1, 1,0 
local, 14, 1,-1, 1,0 
local, 15,0,0,0,0,90,0,0, 

I  BUILD  MODEL 
! CREATE  GEOMETRY 
! ASSIGN  KEYPOINTS 

’NOTE  MODEL  COULD  HAVE  BEEN  BUILT  BY  SYMETRIES 

!  CREATE  MODEL  KEYPOINTS 

csys , 11 

k,  1 , .5,270 

csys , 12 

k,  2 ,  .5,270 

k,  3 ,  .5,0 

csys , 13 

k,4, .5, 0 

k,5, .5, 90 

csys, 14 

k,6, .5, 90 

k, 7, .5,180 

csys , 11 

k ,  8 ,  .5,180 

k, 9, 1,270 

csys , 12 

k, 10, 1,270 

k, 11, 1,0 

csys, 13 

k, 12,1,0 
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k,  13, 1, 90 
csys , 14 
k, 14, 1, 90 
k, 15, 1,180 
csys , 11 

k,  16, 1,180 

! DEFINE  LINES 
csys , 11 

l, 8,1,16 

1.16.9.16 
csys , 12 

1.2.3.16 

1,10,11,16 
csys , 13 

1.4.5.16 

1.12.13.16 
csys , 14 

1.6.7.16 

1.14. 15.16 


csys , 0 

1.1.2.40 

1.9.10.40 

1.3.4.40 

1.11.12.40 

1.5.6.40 

1.13.14.40 

1.7.8.40 

1.15.16.40 

1.1.9.10 

1.2.10.10 

1.3.11.10 

1.4.12.10 

1.5.13.10 

1.6.14.10 

1.7.15.10 

1.8.16.10 

! CREATE  AREAS 

a, 9, 10, 2,1 
a, 10, 11, 3, 2 
a, 11, 12, 4, 3 
a, 12, 13, 5, 4 
a, 13, 14, 6, 5 
a, 14, 15, 7, 6 
a, 15, 16, 8, 7 
a, 16, 9, 1,8 

! ASSIGN  ELEMENT  TYPE,  MATERIALS,  ELEMENT  COORDINATE  SYSTEMS 
! CREATE  MESH 

!  MENU  PATHS  FOR  ESYS : 
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!  Main  Menu  Preprocessor  >Create  >Elements  >Elem  Attributes 
!  or  Main  Menu  Preprocessor  >Define  >Default  Attribs 

esys, 11 
mat ,  1 
type, 1 
amesh, 8 
esys , 12 
mat ,  2 
type , 2 
amesh, 2 
esys, 13 
mat ,  3 
type, 3 
amesh, 4 
esys , 14 
mat ,  4 
type, 4 
amesh, 6 
esys, 15 
mat ,  5 
type, 5 
amesh, 1,5,4 
esys , 0 
mat ,  6 
type, 6 
amesh, 3 ,7,4 

esys ,  0 

!  CONNECT  THE  MESHED  AREA  ELEMENTS  TOGETHER 
!  AND  COMBINE  COMMON  NODES 
nummrg , node 

I  CREATE  (SAVE)  file.db 
save 

!  LEAVE  PREP 7 
f  ini 

! ENTER  SOLUTION  PHASE 
/solu 

'USE  AUTOMATIC  TIME  STEPPING  FOR  NON-LINEAR  PROBLEMS, 
autots, on 

!  CONSTRAIN  STRUCTURE  IN  SPACE 
esys, 0 

nsel , s, loc,y, -1 . 5, 1 . 5 
nsel , r , loc,x, 0 
d, all , ux, 0 
nsel , all 

nsel , s , loc,y, 0 

nsel , r , loc,x, -1 . 5 , 1 . 5 

d, all , uy, 0 
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nsel , all 


!  APPLY  INTERNAL  PRESSURE  TO  STRUCTURE 
csys , 11 

nsel , s, loc,x, . 5 
sf , all , pres , 1000 
nsel , all 

csys , 12 

nsel , s, loc,x, . 5 
sf , all , pres , 1000 
nsel , all 

csys , 13 

nsel , s , loc,x, .5 
sf .all, pres, 1000 
.nsel , all 

csys, 14 

nsel , s , loc , x, . 5 
sf, all, pres, 1000 
nsel , all 

csys , 0 

nsel , s, loc,x,  -1.5, 1 .5 
nsel , r, loc,y, -1,1 
sf, all, pres, 1000 
nsel , all 

csys , 0 

nsel ,s,loc,x,-l,l 
nsel, r, loc,y, -1.5, 1 .5 
sf , all , pres , 1000 
nsel , all 


csys , 0 

1  SOLVE  THE  PROBLEM 
solve 

!  REWRITE  file . db  INCLUDING  SOLUTION  COMMANDS 
save 

!  LEAVE  SOULUTION  PHASE 
f  ini 
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!  ENTER  POST  PROCESSOR 
/postl 

!  INPUT  RESULT  SET 
set ,1,1 

!  PRESENT  RESULTS  IN  ELEMENT  COORDINATE  SYSTEMS 
!  MENU  PATHS  FOR  RSYS : 

!  Main  Menu  >  General  Postproc  >  Options  for  Outp 
1  or  Utility  >  List  >Results>  Options 

rsys, solu 

!  PLOT  AREAS 

aplo 

!  (WAIT  IN  INTERACTIVE  MODE) 

!  BEGIN  POSTPROCESSING  PROCEDURE  FOR  LAMPAT  BY  CREATING 
!  MATID.FIL  AND  STRESS. FIL  FILES. 

!  AFTER  READING  IN  THE  MACRO. WRITE,  POSTPROCESSING  IN  LAMPAT, 

I  AND  WRITING  THE  ETABLE  DATA  WITH  MACRO. READ, 

!  SAVE  THE  DATA  BY  WRITING  TO  file.db  (ENTER  SAVE  IN  THE  COMMAND  LINE 
!  OR  SELECT  THE  SAVE  DB  BUTTON  WITH  THE  MOUSE)  .  THE  ETABLE  DATA  IS 
!  WRITTEN  INTO  THE  file.db,  NOT  THE  file.rst  FILE. 
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Intentionally  left  blank. 


Appendix  E.  Data  Transfer  Macros  (MACRO.WRITE) 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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Intentionally  left  blank. 
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Macro  Write.  This  is  used  to  write  ANSYS  data  into  files  for  use  in  LAMP  AT. 
Associate  materials  and  stresses  with  element  numbers  are  included. 


/com 

/com 

/com  MACRO  FOR  STRESS  ----------- 

/com 

/com 

*get ,  numelem,  elem,  ,  count 

*dim, elenum, array , numelem, 1 , 1 
*do, i , 1 , numelem 
elenum (i) =i 
*enddo 

etable , sx, s ,  x 
etable, sy, s,y 
etable ,sz,s,z 
etable  , syz, s, yz 
etable , sxz , s , xz 
etable , sxy, s, xy 

*dim, esx, array , numelem, 1 , 1 
*dim, esy , array , numelem, 1 , 1 
*dim, esz , array , numelem,  1 , 1 
*dim, esyz, array, numelem, 1, 1 
*dim, esxz , array, numelem, 1, 1 
* dim, esxy, array, numelem,  1, 1 

*vget,esx(l) , elem, 1 , etab, sx 
*vget , esy (1) , elem, 1 , etab, sy 
*vget , esz (1) , elem, 1 , etab, sz 
*vget , esyz (1) , elem, 1 , etab, syz 
*vget , esxz (1) , elem, 1 , etab, sxz 
*vget , esxy (1) , elem, 1 , etab, sxy 

*cfopen, stress , f il 

*vwrite, elenum (1) , esx (1) , esy (1) , esz (1) , esyz (1) , esxz (1) , esxy (1) 

/com  The  following  format  line  may  vary  from  one  LAMPAT  version 

/com  to  another.  Be  sure  to  be  consistent  and  large  enough  for 

/com  your  model . 

(f 6 . 0 , / , 6el3 . 5) 

*cf close 

/com 

/com 

/com  MACRO  FOR  STRAIN  ----------- 

/com 

/com 

/com  Users  might  want  to  remove  the  macro  for 
/com  the  strain  section.  Its  output  was  not 
/com  used  in  this  example.  It  may  consume 
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/com 


a  significant  amount  of  time  for  a  large  model. 


etable , ex, epto, x 
e table, ey, epto, y 
etable , ez , epto, z 
etable , eyz , epto, yz 
etable, exz, epto,xz 
etable, exy, epto,xy 

*dim, epsx, array , numelem, 1 , 1 
*dim, epsy, array, numelem, 1 , 1 
*dim, epsz , array , numelem, 1 , 1 
*dim, epsyz, array, numelem, 1 , 1 
*dim, epsxz , array, numelem, 1 , 1 
*dim, epsxy, array, numelem, 1 , 1 

*vget , epsx (1) , elem, 1 , etab, ex 
*vget , epsy (1) , elem, 1 , etab, ey 
*vget , epsz (1) , elem, 1 , etab, ez 
*vget , epsyz (1) , elem, 1 ,  etab, eyz 
*vget , epsxz ( 1) , elem,  1 ,  etab, exz 
*vget , epsxy (1) , elem,  1,  etab, exy 

*cfopen, strain, f il 

*vwrite, elenum (1) , epsx (1) , epsy  (1) , epsz (1) , epsyz (1) , epsxz (1) , epsxy (1) 
/com  The  following  format  line  may  vary  from  one  LAMPAT  version 

/com  to  another.  Be  sure  to  be  consistent  and  large  enough  for 

/com  your  model . 

(f 6 . 0 , / , 6el3 . 5 ) 

*cf close 

/com 

/com 

/com  MACRO  FOR  MATERIAL  ID  --------- 

/com 

/com 

*cf open,matid, f il 
*do, i , 1 , numelem 
*get ,ml, elem, i , attr,mat 
*vwrite, i , ml 
(2fl3 .5) 

*enddo 
*cf close 


MACRO .  READ 

This  is  used  to  read  LAMPAT  data  into  ANSYS  for  postprocessing.  It  transfers 
stresses,  strains,  failure  mode,  safety  factor,  ply  failure,  and  XXXX  into  element 
data  tables  (etable). 

/com 
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/com 

/com  MACRO  FOR  READING  ------ 

/com 

/com 

*get,numelem,elem, , count 
numt  o  t = nume 1  em*  7 

*dim, elenum, array , nume lem, 1, 1 
*do, i,  1, numelem 
elenum (i) =i 
*enddo 

*dim, a, array, numtot ,1,1 
*dim, al , array, numelem, 1, 1 
*dim, a2 , array, numelem, 1 , 1 
*dim, a3 , array, numelem, 1 , 1 
*dim, a4 , array, numelem, 1 , 1 
*dim, a 5 , array , numelem, 1 , 1 
*dim, a6 , array, numelem, 1 , 1 

*vread,a(l) ,lampat,fil 
/com  The  following  format  line 

/com  to  another.  Be  sure  to  be 

/com  your  model. 

(f 6 . 0 , / , 6el3 .5) 

*do, i , 1 , numelem 
k=  (i-1) *7 
al (i) =a (k+2) 
a2 (i) =a (k+3) 
a3 (i) =a (k+4) 
a4  (i)  =a  (k+5) 
a5 (i) =a (k+6) 
a6 (i) =a (k+7) 

*enddo 

*do, i, 1, numelem 
desol , i , all , eppl , x, al  (i) 
desol, i, all, eppl, y,a2 (i) 
desol , i , all , eppl ,  z ,  a3 ( i) 
desol , i , all , eppl , yz , a4  ( i ) 
desol , i , all , eppl , xz , a5  ( i ) 
desol , i , all , eppl ,xy, a6  (i) 

*enddo 

etable , saf e-f , eppl , x 
e table, f ail-m, eppl , y 
etable, crit-p, eppl , z 
etable , c4 , eppl , yz 
etable , c5 , eppl , xz 
etable , c6 , eppl , xy 


may  vary  from  one  LAMPAT  version 
consistent  and  large  enough  for 
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Macro.read  uses  the  standard  LAMPAT  code  processing  data  file,  which 
is  named  lampat.fil  prior  to  executing  this  macro.  This  macro  writes  information 
(safety  factor  data,  failure  mode  data,  and  ply  failure  data)  into  etable  data  as 
plastic  strain.  To  examine  the  plastic  strain  in  the  model,  this  must  be  done  prior 
to  executing  the  macro.  The  following  is  an  alternative  macro  that  uses  the 
ANSYS  smisc  data  terms  as  portals  for  the  LAMPAT  results. 


/com 

/com 

/com  MACRO  FOR  READING  -  -  - 

/com 

/com 

*get , numelem, elem, , count 
numt  o t = nume 1 em*  7 

*dim, elenum, array , numelem, 1 , 1 
*do, i; 1, numelem 
elenum (i) =i 
*enddo 

*dim, a , array , numt ot ,1,1 
*dim, al , array, numelem, 1 , 1 
*dim, a2 , array , numelem, 1 , 1 
*dim, a3 , array, numelem, 1 , 1 
*dim, a4 , array, numelem, 1 , 1 
*dim, a5 , array, numelem, 1 , 1 
*dim, a6 , array , numelem, 1 , 1 


*vread , a ( 1 ) , lampat . f i 1 

/com  The  following  format  line  may  vary  from  one  LAMPAT  version 

/com  to  another.  Be  sure  to  be  consistent  and  large  enough  for 

/com  your  model . 

(f 6 . 0 , / , 6el3 . 5) 

*do, i , 1 , numelem 
k=  (i-1) *7 
al (i) =a (k+2 ) 
a2 (i) =a (k+3 ) 
a3 (i) =a (k+4 ) 

/com  Since  a4,  a5,  and  a 6  are  currently  not  used,  they  are  not  read, 

/com  To  save  time,  remove  them  in  the  LAMPAT  code  to  save  time 

/com  when  running  large  models. 

*enddo 

/com  The  selection  of  the  variable  term  used  to 
/com  read  data  into  the  element  data  tables  in 
/com  ANSYS  varies  with  the  element  type  used  (i.e. 

/com  SOLID  42) .  If  the  user  wants  to  get 

/com  results  on  the  smisc  data,  other  terms  should 

/com  be  used.  Alternatively,  if  the  user  postprocesses 
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/com  the  ANSYS  data  (i.e.  smisc  values) 

/com  prior  to  reading  in  the  lampat  results  (that 
/com  is  prior  to  executing  macro. read),  the  results 

/com  will  not  be  corrupted  by  this  read  in  data. 

etable, safe-f , smisc, 1 
etable, fail-m, smisc, 2 
e table, crit-p, smisc , 3 

*do, i, l,numelem 
detab, i, safe-f , al (i) 
detab, i, fail-m, a2 (i) 
detab, i, crit-p, a3 (i) 

*enddo 
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Intentionally  left  blank. 


NO.  OF  NO.  OF 

COPIES  ORGANIZATION  COPIES  ORGANIZATION 


V 


2  DEFENSE  TECHNICAL 

INFORMATION  CENTER 
DTIC  OCA 

8725  JOHN  J  KINGMAN  RD 
STE  0944 

FT  BELVOIR  VA  22060-6218 

1  HQDA 

DAMO  FDT 

400  ARMY  PENTAGON 

WASHINGTON  DC  20310-0460 


3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  Cl  LL 
2800  POWDER  MILL  RD 
ADELPHI MD  20783-1197 

3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  Cl  1ST 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 


1  OSD 

OUSD(A&T) /  ODDR&E(R)  ABERDEEN  PROVING  GROUND 

DRRJTREW 

3800  DEFENSE  PENTAGON  2  DIR  USARL 

WASHINGTON  DC  20301-3800  AMSRL  Cl  LP  (BLDG  305) 


1  COMMANDING  GENERAL 

US  ARMY  MATERIEL  CMD 
AMCRDATF 
5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333-0001 


1  INST  FOR  ADVNCD  TCHNLGY 

THE  UNIV  OF  TEXAS  AT  AUSTIN 
3925  W  BRAKER  LN  STE  400 
AUSTIN  TX  78759-5316 


1  US  MILITARY  ACADEMY 

MATH  SCI  CTR  EXCELLENCE 
MADN  MATH 
THAYER  HALL 
WEST  POINT  NY  10996-1786 


1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  D 
DR  D  SMITH 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  Cl  AI R 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 
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NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  CP  CA 
D  SNIDER 

2800  POWDER  MILL  RD 
ADELPHI MD  20783-1145 

1  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  Cl  IS  R 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 

3  DIRECTOR 

US  ARMY  RESEARCH  LAB 
AMSRL  OP  SD  TL 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 

1  DPTY  ASST  SECY  FOR  R&T 
SARD  TT 
THE  PENTAGON 
RM  3EA79 

WASHINGTON  DC  20301-7100 

1  COMMANDER 

US  ARMY  MATERIEL  CMD 
AMXMI INT 

5001  EISENHOWER  AVE 
ALEXANDRIA  VA  22333-0001 

4  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  CC 

G  PAYNE 
J  GEHBAUER 
C  BAULIEU 
HOPAT 

PICATINNY  ARSENAL  NJ 
07806-5000 

2  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  AE  WW 
E  BAKER 

J  PEARSON 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  FSE 
PICATINNY  ARSENAL  NJ 
07806-5000 


1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  AR  TD 
C  SPINELLI 

PICATINNY  ARSENAL  NJ 
07806-5000 

6  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  CCH  A 
W  ANDREWS 
S  MUSALLI 
R  CARR 
M  LUCIANO 
E  LOGSDEN 
T  LOUZEIRO 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  CCH  P 
J  LUTZ 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  AR  FSF  T 
C  LIVECCHIA 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 

US  ARMY  ARDEC 
AMSTA  ASF 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  QACTC 
C  PATEL 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 
US  ARMY  ARDEC 
AMSTA  ARM 
D  DEMELLA 

PICATINNY  ARSENAL  NJ 
07806-5000 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


2  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  FSP  G 
M  SCHIKSNIS 
D  CARLUCCI 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  FSP  A 
P  KISATSKY 

PICATINNY  ARSENAL  NJ 
07806-5000 

2  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  CCH  C 
H  CHANIN 

S  CHICO 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  QACT 
D  RIGOGLIOSO 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  WET 
TSACHAR 
BLDG  172 

PICATINNY  ARSENAL  NJ 
07806-5000 

1  US  ARMY  ARDEC 

INTELLIGENCE  SPECIALIST 
AMSTA  AR  WELF 
M  GUERRIERE 
PICATINNY  ARSENAL  NJ 
07806-5000 


9  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  CCH  B 
P  DONADIA 
F  DONLON 
P  VALENTI 
C  KNUTSON 
G  EUSTICE 
S  PATEL 
G  WAGNECZ 
RSAYER 
F  CHANG 

PICATINNY  ARSENAL  NJ 
07806-5000 

6  COMMANDER 
US  ARMY  ARDEC 
AMSTA  AR  CCL 
F  PUZYCKI 
R  MCHUGH 
D  CONWAY 
E  J  AROSZEWSKI 
R  SCHLENNER 
M  CLUNE 

PICATINNY  ARSENAL  NJ 
07806-5000 

5  PM  ARMS 

SFAEGCSSSD 
COL  SLEDGE 
M  DEVINE 
W  DEMASSI 
J  PRITCHARD 
S  HROWNAK 
PICATINNY  ARSENAL  NJ 
07806-5000 

2  PEO  FIELD  ARTILLERY  SYS 

SFAE  FAS  PM 
H  GOLDMAN 
T  MCWILLIAMS 
PICATINNY  ARSENAL  NJ 
07806-5000 

1  COMMANDER 
US  ARMY  ARDEC 
AMSTA  ARWEA 
J  BRESCIA 

PICATINNY  ARSENAL  NJ 
07806-5000 


NO.  OF 

COPIES  ORGANIZATION 

1  DIRECTOR 

AIR  FORCE  RESEARCH  LAB 
MLLMD 
D  MIRACLE 
2230  TENTH  ST 

WRIGHT  PATTERSON  AFB  OH 
45433-7817 

2  COMMANDER 

US  ARMY  AMCOM 
AVIATION  APPLIED  TECH  DIR 
J  SCHUCK 

FT  EUSTIS  VA  23604 -5577 

1  BENET  LABORATORIES 
AMSTA  AR  CCB 
EKATHE 

WATERVLIETNY  12189-4050 

1  DIRECTOR 

US  ARMY  AMCOM 
SFAE  AV  RAM  TV 
D  CALDWELL 
BLDG  5300 

REDSTONE  ARSENAL  AL 
35898 

1  NAVAL  SURFACE  WARFARE  CTR 
DAHLGREN  DIV  CODE  G06 
DAHLGREN  VA  22448 

1  US  ARMY  COLD  REGIONS 
RSCH  &  ENGRNG  LAB 
PDUTTA 
72  LYME  RD 
HANOVER  NH  03755 

1  USA  SBCCOM  PM  SOLDIER  SPT 
AMSSB  PM  RSS  A 

J  CONNORS 
KANSAS  ST 

NATICK  MA  01760-5057 

2  OFC  OF  NAVAL  RESEARCH 
D  SIEGEL  CODE  351 

J  KELLY 

800  N  QUINCY  ST 
ARLINGTON  VA  22217-5660 


NO.  OF 

COPIES  ORGANIZATION 

1  NAVAL  SURFACE  WARFARE  CTR 
TECH  LIBRARY  CODE  323 
17320  DAHLGREN  RD 
DAHLGREN  VA  22448 

1  NAVAL  SURFACE  WARFARE  CTR 
CRANE  DIVISION 

M  JOHNSON  CODE  20H4 
LOUISVILLE  KY  40214-5245 

2  NAVAL  SURFACE  WARFARE  CTR 
U  SORATHIA 

C  WILLIAMS  CD  6551 
9500  MACARTHUR  BLVD 
WEST  BETHESDA  MD  20817 

2  COMMANDER 

NAVAL  SURFACE  WARFARE  CTR 
CARDEROCK  DIVISION 
R  PETERSON  CODE  2020 
M  CRITCHFIELD  CODE  1730 
BETHESDA  MD  20084 

1  NAVAL  SEA  SYSTEMS  CMD 
D  LIESE 

2531  JEFFERSON  DAVIS  HWY 
ARLINGTON  VA  22242-5160 

1  NAVAL  SURFACE  WARFARE  CTR 
M  LACY  CODE  B02 
17320  DAHLGREN  RD 
DAHLGREN  VA  22448 

1  EXPEDITIONARY  WARFARE 

DIV  N85 
F  SHOUP 

2000  NAVY  PENTAGON 
WASHINGTON  DC  20350-2000 

8  NAVAL  SURFACE  WARFARE  CTR 

J  FRANCIS  CODE  G30 
D  WILSON  CODE  G32 
R  D  COOPER  CODE  G32 
J  FRAYSSE  CODE  G33 
E  ROWE  CODE  G33 
T  DURAN  CODE  G33 
L  DE  SIMONE  CODE  G33 
R  HUBBARD  CODE  G33 
DAHLGREN  VA  22448 
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NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


t 


2  NAVAL  SURFACE  WARFARE  CTR 
CARDEROCK  DIVISION 
R  CRANE  CODE  2802 
C  WILLIAMS  CODE  6553 
3A  LEGGETT  CIR 
BETHESDA  MD  20054-5000 

1  AFRL  MLBC 

2941  P  ST  RM  136 

WRIGHT  PATTERSON  AFB  OH 

45433-7750 

1  AFRL  MLSS 
R THOMSON 

2179  12TH  ST  RM  122 
WRIGHT  PATTERSON  AFB  OH 
45433-7718 

2  AFRL 

F  ABRAMS 
J  BROWN 
BLDG  653 
2977  P  ST  STE  6 

WRIGHT  PATTERSON  AFB  OH 
45433-7739 

1  WATERWAYS  EXPERIMENT 
D  SCOTT 

3909  HALLS  FERRY  RD  SC  C 
VICKSBURG  MS  39180 

1  OSD 
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